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HO 要 :针对 沥青 质 颗 粒 在 水 平 并 并 简 的 沉积 问题 ,本 研究 采用 Fluent 中 的 DPM 模块 追踪 沥青 质 
颗粒 在 油管 中 的 运动 轨迹 ,采用 realizable k-e 模型 对 尺寸 为 10 ~ 100 um 的 沥青 质 颗粒 在 水 平井 简 
内 壁 的 沉积 进行 了 研究 。 分 析 了 重力 、Saffman 力 以 及 阻力 作用 下 , 稀 度 、 颗 粒 尺 寸 以 及 流速 对 沥青 
质 颗粒 沉积 速率 的 影响 。CFD 模拟 结果 表明 : pli 沉积 速率 随 着 黏度 的 增 大 而 减 小 ,这 主 
要 是 因为 层 流 阶段 ,流体 流动 较 平 稳 , 使 得 管 壁 附近 速度 梯度 减 小 ,Safftman 升力 降低 引起 的 ， in 
3 LUE, 3: 45A 0.0012 Pa * s 增加 至 0.012 Pa * s 4 沥青 质 沉积 速率 会 缓慢 增加 ; 4A 
0.012 Pa * s 增加 至 0.053 Pa * s, 沉 积 速 率 和 沉积 百分数 会 随 之 下 降 ;@) 当 颗粒 直径 为 10 ~ 100 um 
时 ,沥青 质 颗粒 沉积 百分数 和 沉积 速率 会 随 着 流速 的 增加 而 减 小 ;加 当 速度 为 0.3 ~2 m/s 时 ,沥青 
质 颗粒 沉积 百 分 E rd a 研究 结果 为 水 平 并 井 简 中 的 沥青 沉 


行为 指明 了 方向 ,认为 合理 控制 并 简 中 流体 速度 、 黏 度 和 颗粒 直径 可 作为 去 除 并 简 沥 青 质 沉积 的 有 
效 手 段 。 
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Numerical simulation of asphaltene particle transport in 
horizontal well horizontal section based on DPM model 
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Abstract :To study the problem of asphaltic particles deposition in horizontal wellbore ,the DPM module in 
Fluent was used to track the particle movement in tubing. And the realizable k-e model was adopted to in- 
vestigate the asphaltene aggregates with a size of 10 ~ 100 jum on the inner wall of tubing. The effects of 
particle size and velocity on the asphaltene particle deposition rate are analyzed when the gravity , saffman 
lift, and drag forces are exerted on the particles. CFD simulation results show: (I) The deposition rate de- 


creases with the increase of viscosity in the laminar flow stage. This is mainly because in the laminar flow 
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stage ,the fluid flow is relatively stable and Saffman lift decreases. The asphaltene deposition rate increases 


slowly in the range from 0. 001 2 Pa - 


the viscosity increases from 0.012 Pa - 


s to 0.012 Pa * 
s to 0.053 Pa - 


s;the deposition rate and percentage decrease when 


s in the turbulent stage. (2)The deposition percent- 


age and deposition rate of asphaltene particles decrease with the increase of flow rate when the particle di- 


ameter is 10-100 jum. 3) When the velocity is 0. 3-2 m/s ,the percentage and deposition rate of asphaltene 


particles increase with the increase of particle diameter. The analysis shows that reasonable control of well- 


bore fluid velocity and particle diameter can be an effective means to reduce asphaltic deposition in well- 


bore. 
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在 石油 工业 生产 过 程 中 ,由 于 原油 温度 .压力 以 
及 组 分 的 改变 ,溶解 在 油 中 的 沥青 质 稳定 性 会 被 打 
破 ,沥青 质 沉淀 会 从 油 流 中 析出 ,部 分 沥青 质 颗粒 会 
黏附 在 井 简 油管 内 壁 ,堵塞 油 流 的 有 效 通道 ,严重 降 
候 肖 井 产能 , 极 大 地 增加 了 油田 的 开采 成 本 A 
此 研究 井 简 沥 青 质 颗粒 的 沉积 特性 变 得 至 关 重要 。 
ST ESCOBEDO 等 1 在 建立 固体 颗粒 沉积 模型 时 ， 
考 感 了 布朗 亲 流 扩散 系数 以 及 沥青 质 颗粒 输送 到 
管 忌 的 惯性 效应 ,以 此 预测 沥青 质 颗 粒 在 清流 条 件 
KIRE, RAMIREZ 等 “基于 分 子 扩散 和 剪 
切 海 沥青 质 颗粒 径 向 扩散 竞争 机 制 的 假设 ,利用 多 
相 阁 组 分 流体 动力 学 模型 预测 生产 井 的 沥青 质 颗 粒 
沈 租 ,该 模型 预测 了 沥青 质 沉积 质量 剖面 随时 间 和 
位 置 的 变化 规律 ,同时 还 认为 径 向 温度 梯度 是 沥青 
质 颗粒 运 移 到 管 壁 的 原因 。 但 是 ESKIN'"' , JAMI- 
ALAHMADI ? 等 认为 沥青 质 颗粒 在 管 壁 上 的 沉积 速 
率 是 由 化 学 反应 机 理 控制 ,而 并 非 是 温度 梯度 的 控 
制 5 用 在 方程 中 加 入 指数 项 以 考虑 温度 对 化 学 反应 
figu, "JE tg og Ut 认为 沥青 质 颗 粒 向 井 壁 运 移 
的 动力 主要 是 分 子 扩散 力 和 惯性 力 ,并 且 研 究 了 流 
速 , 粒 径 等 参数 对 沉积 速率 的 影响 。SOULGANIL 
建立 了 反应 ,扩散 受 限 聚集 的 沥青 质 聚 集 动力 学 的 
通 式 ,研究 了 重力 .阻力 、 浮 力 对 沥青 质 颗粒 平衡 作 
JH, SEYYEDBAGHERI 4?! fi Hj Fluent 软件 中 Eu- 
lerian-Lagrangian 方法 模拟 了 3 维 油 管 中 沥 青 质 聚集 
颗粒 的 沉积 特性 ,模拟 结果 显示 ,沥青 质 颗 粒 沉积 速 
度 与 颗粒 数目 以 及 管 壁 粗 糙 度 无 关 , 沉积 通 量 会 随 
着 流速 的 增加 而 降低 。 综 上 所 述 , 以 往 的 数值 研究 
大 多 针对 的 是 垂直 井 简 内 管 壁 沥 青 质 聚集 体 的 沉积 
行为 , 鲜 少 有 关于 水 平井 井 简 的 颗粒 沉积 行为 的 研 
究 报道 ;另外 学 者 们 在 用 CFD 软件 模拟 过 程 沥青 质 
沉积 过 程 时 ,并 未 考虑 原油 黏度 对 沥青 质 沉 积 速率 
的 影响 。 


为 此 ,本 研究 采用 Lagranian 方法 建立 沥青 质 聚 
集体 的 沉积 数学 模型 ,利用 单 向 耦合 的 方法 研究 水 
平井 井 简 内 壁 的 沥青 质 颗粒 沉积 ,其 中 采用 DPM 模 
块 追 踪 油 管 中 颗 粒 运动 轨 迹 , 利 用 realizable k-e 模 
型 模拟 井 简 中 充分 发 展 的 灌流 条 件 。 研 究 了 颗粒 直 
径 .速度 以 及 黏度 对 沥青 质 沉 积 百分数 以 及 沉积 速 
率 的 影响 ,该 研究 结果 对 水 平井 井 简 沥青 质 沉积 的 
去 除 具 有 指导 意义 。 


1 数值 模型 描述 


1.1 沥青 质 聚 集体 沉积 数学 模型 


依据 Lagranian 方法 跟踪 颗粒 呈 ] ,模拟 沥青 质 
颗粒 在 水 平井 简 中 的 运动 情况 ,对 于 沥青 质 颗粒 的 
性 质 作 以 下 假设 :四 沉淀 运动 过 程 中 ,密度 保持 不 
变 ;@ 颗 粒 为 刚体 ,不 会 发 生变 形 ;@ 颗 粒 为 质点 , 忽 
略 颗粒 -颗粒 之 间 的 相互 作用 。 根 据 Eulerian-La- 
grangian 模型 ,建立 液 相 的 连续 方程 与 动量 守恒 
Vir 。 

油 相 连续 性 方程 为 


Gp.) + V: (aup) 20 (1) 
油 相 动量 方程 为 
全 (auuup。) 二 
-a, Vp + Vv, tapB th (2) 


式 中 :a 为 液 相 的 体积 浓度 分 数 ,% ;p, 表 示 油 相 的 
密度 ;V 表 示 哈 密 顿 算 子 ;wu, 表 示 油 相 的 速度 , m/s;p 
表示 压强 ,Pa;7, 表示 应 力 张 量 ; 了 表示 两 相间 的 体 
积 作 用 力 ,N。 

两 相间 的 体积 作用 力主 要 包含 重力 浮力. 阻 
力 、Saffman 升力 、 虚 拟 质 量力 、 压 力 梯 度 力 、 静 电力 、 
Basset 力 和 Magnus 力 等 ,由 于 后 面 5 个 作用 力 在 数 
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量 级 上 相对 较 小 ,一 般 可 忽略 。 颗 粒 使 用 Lagra- 
nian 轨道 模型 ,但 当 颗 粒 浓度 稀薄 时 (体积 分 数 < 
10% ~12% ) ,只 考虑 流 场 对 颗粒 的 作用 。 因 此 , 根 
据 牛顿 第 二 定律 ,单个 颗粒 运动 的 动力 学 沉积 
程 为 


m $ 

P dt 

式 中 :Fi 表示 沥青 质 颗 粒 的 阻力 ,N;G 表示 重力 ,N; 

mm 表示 沥青 质 聚 集体 的 质量 ,kg;v 表示 流体 的 速 

JE m/s; F, zn. Saffman 升力 , N; 下 标 p 表示 单个 

颗粒 。 

1) 阻 力 

沥青 质 颗粒 在 流动 过 程 中 受到 的 阻力 是 由 颗粒 

表面 的 压力 和 剪 切 力 引起 的 ”。 阻力 通常 与 颗粒 

斥 污 形状 以 及 流速 大 小 及 方向 相关 ,其 相关 式 可 以 


=G+F,+F (3) 


示 为 

Op Lp Gi - uj)" Cs (4) 
= (70 8G 
Xr C, 4i Cunningham f E cs d, t f 


粹 局 寸 ,m;w 指 的 是 7 方向 流体 的 速度 ,m/s;u, 表 示 
r 力 向 颗粒 的 速度 ,m/s;Cb 表 示 自力 系数 。 


ARIRE URRH 
N =, Re, <1 
Pra 5) 


a2) 24 (1 +0. 158,9), 1<Re,<1000 ( 
* | Re p I 


© F 
E 0. 44 Re, > 1 000 
E2) Saffman 升力 

昌 油 流 中 的 沥青 质 颗粒 由 于 剪 切 作用 ,分 散 的 沥 
青 质 颗 粒 会 受到 垂直 于 流动 方向 的 升力 ,促使 颗粒 
获得 动量 , 随 着 油 一 起 运动 "” 。 利 用 Saffman 方程 
计算 了 剪 切 力 对 颗粒 的 作用 大 小 。 

F, = 1.615p "di (y - 9) |y |? (6) 

其 中 :v 为 流体 的 运动 速度 ,m/s;y 表示 沥青 质 颗 粒 
在 随 油 流动 过 程 中 所 受到 的 剪 切 速率 ,m/s。 


1.2 物理 模型 


本 小 节 人 研究 的 水 平井 简 模 型 来 自 文献 [20 ] ,内 
径 70 mm ,长 度 2 000 mm, 具 体 模 型 如 图 1 所 示 。 由 
于 原油 流动 是 在 充分 发 展 的 油 流 条 件 下 进行 的 ,为 
消除 入 口 效 应 的 影响 ,首先 在 足够 长 的 三 维 水 平井 
简 中 模拟 原油 流动 ,然后 将 输出 结果 用 于 人 研究 水 平 
井 简 的 入 口 油 流 的 初始 阶段 。 
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由 于 本 研究 主要 是 模拟 沥青 质 颗 粒 在 油管 壁面 
的 沉积 情况 ,因此 为 了 更 真实 地 模拟 管 壁 的 复杂 流 
场 以 及 沥青 质 颗粒 在 管 壁 处 的 沉积 过 程 ,本 研究 利 
用 Proximity and Curvature 方法 在 油管 壁面 设置 了 
10 层 精 细 网 格 ,节点 数 约 72 万 ,网 格 数 约 为 182 H, 


既 能 保证 连续 油 流 相 的 仿真 精度 ,又 可 避免 因 网 格 
尺度 偏 小 引起 DPM 模型 计算 出 错 的 问题 ,三 维 水 平 
井 简 物 理 模型 如 图 2 所 示 。 
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图 1 物理 模型 
Fig.1 Physical model 


(a) 主 视图 


到 2 水 平井 简 物 理 模型 
Fig.2 The physical model of horizontal wellbore 

为 了 对 网 格 数量 进行 独立 性 验证 ,本 研究 将 网 
格 数量 增加 至 511 万 ,统计 了 油管 内 壁 沥青 质 颗粒 
的 沉积 数量 ,将 其 沉积 数目 与 182 万 网 格 时 的 沥青 
质 颗粒 进行 对 比 ,如 表 1 所 示 。 由 表格 可 看 出 ,相对 
于 182 万 网 格 , 误 差 稳定 在 6.1% 左右 ,考虑 到 建 模 
时 ,运行 时 间 随 网 格 尺寸 的 增加 而 显著 增加 , 故 本 研 
究 选 择 节 点 数 为 72 万 ,网 格 数 为 182 万 的 模型 。 

表 1 网 格 数目 敏感 性 分 析 


Tab.1 The sensitivity analysis of mesh number 


(b) 入 口 截面 


节点 数 网 格 数 沉积 颗粒 数 误差 /% 
268 064 758 342 5717 36.6 
913 259 3 074 764 4 448 6.1 

1 430 233 5115 943 4041 6.2 


1.3 边界 条 件 设 置 


本 研究 在 计算 过 程 中 , 消 流 模型 采用 Realizable 
k-e 模型 ,固体 颗粒 的 追踪 采用 Fluent 中 的 DPM 离 
散 相模 块 ,沥青 质 聚 集 颗 粒 沿 人 口 垂 直面 均匀 注入 ， 


速度 与 油 流 的 速度 保持 一 致 。 在 模拟 过 程 中 主要 考 
虑 重力 .阻力 以 及 Saffman 升力 ,入 口 出口 设 为 逃逸 
型 边界 (escape) , 管 壁 设 为 捕 集 型 边界 (trap ) 。 另 
外 ,人口 压力 为 30 MPa ,出 口 压力 为 29.5 MPa ;温度 
为 50 CC 。 其 中 , 油 为 连续 相 , 沥 青 质 颗粒 为 分 散 相 ， 
详细 参数 如 表 2 所 示 。 

R2 模拟 参数 设 定 


Tab.2 The simulation parameters 


参数 数值 
原油 密度 / (kg * m73) 866 
沥青 质 密度 / (kg : m^?) 1 100 
原油 黏度 /(Pa . s) 0.0012,0.012 ,0.06,0.12,1.2 


入 口 流速 /(m :ss-!) 
颗粒 直径 /um 


0.3,0.5,1.0,1.5,2.0 


10 ~ 100 


1) 颗粒 沉积 速率 


<J 颗 粒 沉积 速率 ,是 描述 固体 表面 颗粒 沉积 特性 
效 物 理 参数 ,其 表达 式 为 “” 
z J _ UnemA N inject 
n= DL uA ) (7) 


JE: /表示 为 单位 时 间 内 颗粒 运 移 到 管 壁 的 质量 
EH C. 表示 颗粒 的 浓度 ,kg/m; PU NEZ | 流 的 
PEE m/s; ;4 表示 井筒 的 横 截 面积 ,m ; 表示 
APRI PES, m; L RRHH JE ms Nu 表示 水 平 
JER EAE D Y TCR LAO I ;Nai 表示 人 
口 段 注 入 的 沥青 质 聚 集体 的 个 数 。 

通常 用 无 量 纲 的 颗粒 沉积 速度 和 无 量 纲 并 驳 时 
间 标 准 颗 粒 的 沉积 速度 '” ,其 中 无 量 纲 表达 式 为 

J 
"PE (8) 
其 中 ut 表示 流体 剪 切 速度 ,可 由 式 (9) 来 计算 。 
w= Unn aE (9) 

其 中 /表示 范 宁 常数 ,可 以 表示 为 


二 =-3.6log[2 + ( k, ) | (10) 


过 
V, z 


J Re 3.7D 
沥青 质 颗粒 沉积 通 量 可 用 下 式 来 计算 。 
J = VS,C, (11) 
2 ) 颗粒 沉积 百分数 
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颗粒 沉积 百分数 表示 水 平井 简 管 壁 上 沥青 质 案 
集 颗 粒 捕获 情况 。 其 表达 式 为 


Nas 
n = Ut x 100% (12) 


另外 ,本 研究 Fluent 模拟 过 程 中 ， 认为 残 差 曲线 
控制 住 10 一 以 下 时 ,为 合理 的 范围 内 ,属于 收敛 。 


2.1 流 场 分 析 


1) 速度 场 分 布 

本 研究 通过 数值 模拟 分 析 了 流速 为 2.0 m/s, 颗 
粒 直径 范围 为 10 ~ 100 pm 时 , 油 流 在 直径 为 
70 mm ,长 为 2 000 mm 的 水 平井 简 内 速度 分 布 情况 ， 
其 结果 如 图 3 所 示 。 


Vector-2 

Velocity Magnitude 
6.22e-01 8.76e-01 

[m/s] 


1.13e+00 1.38e+00 1.64e+00 1.89e+00 2.14e+002.31e+00 


图 3 水 平井 简 内 速度 场 分 布 
Fig.3 Schematic of velocity field distribution 

从 图 3 可 以 看 出 ,水 平井 简 中 的 油 样 在 径 向 流 
动 过 程 中 ,流体 速度 是 呈正 态 分 布 趋 势 的 ,具体 表现 
为 井 简 中 间 油 流速 度 最 大 ,壁面 两 端的 速度 最 小 , 且 
随 着 流体 入 口 的 速度 增加 , 井 简 内 的 最 大 流速 会 向 
管 壁 延伸 ,致使 管 壁 附近 出 现 较 大 的 速度 梯度 , 越 靠 
近 出 口 速度 梯度 越 明 显 。 这 是 因为 随 着 速度 的 增 
加 ,颗粒 在 管 壁 附近 受到 Saffman 升力 会 增加 ,从 而 
会 出 现 速度 梯度 增 大 。 

2) 压 力 场 分 布 

图 4 为 油 样 流速 为 0.5 m/s, 沥 青 质 颗粒 尺寸 为 
10 ~ 100 pm 时 水 平井 简 内 部 的 压力 分 布 图 。 从 图 
中 可 以 看 出 ,沥青 质 颗粒 随 油 流 从 入 口 运 动 至 出 口 
时 ,整个 水 平井 简 的 压力 下 降 速 度 变 化 很 小 ,这 表明 
压力 对 流 场 的 影响 作用 较 小 。 


contour-2 
Static Pressure 
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图 4 水 平井 简 内 压力 场 分 布 
Fig.4 Schematic of pressure field distribution 

3) 颗粒 轨迹 

图 5 为 流速 为 2.0 m/s ,颗粒 直径 范围 为 10 ~ 
100 jum 时 ,沥青 质 颗粒 在 水 平井 简 内 的 颗粒 追踪 轨 
迹 图 。 从 图 中 可 以 看 出 ,由 于 灌流 作用 ,一 部 分 壁面 
的 颗粒 会 在 壁面 发 生 反 弹 ,从 而 偏离 原来 的 运动 轨 
迹 ; 另 一 部 分 沥青 质 颗粒 会 沉积 到 水 平井 简 壁 面 ( 红 
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色 方 框 标 记 ) ,经 过 一 段 时 间 ,沉积 层 厚度 会 逐渐 累 
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积 ,最 终 部 分 堵塞 油 流通 道 。 


6.76e-01 8.45e-01 1.01e*00 1.13e-00 


到 5 水 平井 简 内 颗粒 运动 分 布 


Fig.5 Schematic of particle moving field distribution 


2.2 黏度 对 颗粒 沉积 速率 的 影响 


夭 度 是 通过 雷诺 数 的 变化 来 影响 井 简 中 沥青 质 
颗粒 沉积 的 。 本 研究 分 别 模拟 了 流速 为 0. 5 m/s, W 
粒 尺 寸 为 10 ~ 100 um 时 ,5 PPA 0.001 2,0.012, 
0306 0. 12.1.2 Pa - 。 下 沥青 质 颗粒 沉积 ,其 沉积 
分 亚 和 沉积 速率 如 图 6 ~ 图 7 所 示 。 
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ig.6 The effect of velocity on deposition percentage 
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到 7 黏度 对 沥青 质 颗粒 沉积 速率 的 影响 
Fig.7 The effect of velocity on deposition rate 

在 层 流 阶段 (w => 0.053 Pa . s, Rex2 300) , DC 
积 百 分 数 和 速率 随 着 黏度 的 增 大 而 减 小 ,这 是 因为 
随 着 原油 黏度 的 增加 ,流体 的 流动 逐渐 变 得 平稳 , 管 
壁 附近 的 流体 梯度 会 逐渐 降低 ,致使 Saftman 升力 逐 
渐 减 小 ,从 而 导致 沉积 速率 降低 。 在 潮流 阶段 (Re > 
2 300) ,原油 黏度 从 0. 001 2 Pa.s 增 大 至 0. 012 Pa.s 
时 ,油管 中 原油 状态 为 强 油 流 ,沉积 百分数 和 沉积 速 
率 会 随 之 增 大 , 这 是 因为 原油 黏度 与 油 样 中 的 沥青 


e 
- 


e 
io 


沉积 速率 /(10” *$ In ; 
o 
oo 


质 含量 呈 递 减 指数 关系 , 黏度 越 小 ,原油 中 沥青 
质 含量 越 少 ,从 而 致使 沥青 质 沉 淀 量 减 小 。 当 原油 
黏度 从 0.012 Pa - s 降 至 0.053 Pa * s 时 ,Re 会 不 断 
减 小 (从 10 103 降 至 2300) ,致使 油管 中 的 原油 的 惯 
性 力 减 小 ,使 得 沉积 速率 和 沉积 百分数 降低 。 


2.3 颗粒 尺寸 对 沉积 速率 的 影响 


在 流速 为 2 mys ,黏度 为 0.0012 Pa * s 时 ,本 研 
究 分 别 模拟 沥青 质 颗 粒 直径 为 10、30、50、70、 
100 pm 时 水 平井 简 内 颗粒 运 移 情况 ,从 而 得 到 不 同 
颗粒 尺寸 下 的 沉积 情况 ,如 图 8 所 示 。 从 图 8 可 以 
看 出 , 随 着 沥青 质 颗 粒 尺 寸 的 增加 ,沉积 百分数 随 之 
增加 ,这 是 因为 颗粒 粒 径 的 增 大 ,使 得 颗粒 的 重力 以 
及 Saffman 升力 越 来 越 大 ,颗粒 在 水 平井 简 壁 面 沉 积 
的 趋势 会 大 大 增加 。 图 9 为 5 种 颗粒 尺寸 的 沉积 速 
率 图 ,从 图 中 可 以 看 出 , 当 颗 粒 直径 越 大 ,沉积 速率 
越 大 ,这 表明 重力 是 大 粒 径 颗粒 的 主导 力 。 


0.54 


0s 20 40 60 80 100 120 
颗粒 尺寸 /pm 
图 8 颗粒 尺寸 对 沉积 百分数 的 影响 
Fig.8 The effect of particle size on deposition percentage 
为 了 进一步 验证 上 述 结 论 ,本 研究 在 水 平井 人 
口 端 注 入 直径 范围 为 10 ~ 100 pum 沥青 质 颗粒 ,通过 
分 析出 口 颗粒 的 直径 分 布 ,从 而 确定 哪 种 颗粒 沉积 
的 数量 最 多 。 图 10 为 出 口 颗粒 尺寸 分 布 的 饼 图 。 
从 图 10 中 可 以 看 出 ,出 口 沥青 质 颗粒 直径 范围 主要 
为 10 ~80 pm 的 颗粒 ,这 表明 90 ~ 100 hm 的 沥青 质 
颗粒 大 部 分 都 沉积 在 了 油管 壁面 。 该 结论 与 张 丽 珍 
等 "研究 气 固 两 相 流 时 得 出 的 结论 一 致 , 即 沥青 质 


颗粒 直径 越 大 ,沉积 越 严重 。 
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Fig.9 The effect of particle size on deposition rate 
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Fig. 10 The particle size distribution in the outlet 
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2 流速 对 颗粒 沉积 速率 的 影响 


_ 伍 本 研究 分 别 模拟 了 流速 为 0.3、0.5、1.0、1.5、 
220 m/s 时 ,直径 范围 为 10 ~ 100 um 的 沥青 质 颗粒 
ARFA H HEM, WE 11 ~ 图 12 所 示 。 
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图 11 流速 对 沉积 百分数 的 影响 
Fig.11 The effect of oil velocity on deposition percentage 
从 图 可 以 看 出 ,沥青 质 颗粒 沉积 百分数 和 沉积 

速率 会 随 着 流速 的 增加 而 减 小 。 这 是 因为 随 着 油 速 
的 增加 ,阻力 也 会 随 之 增加 ,壁面 会 产生 较 大 的 前 切 
应 力 , 从 而 导致 沥青 质 颗粒 在 壁面 的 恭 附 概率 降低 ， 
同时 也 会 剥离 管 壁 上 已 沉积 的 沥青 质 。 该 结论 与 
ESKIN 等 ”用 Real View 研究 计算 出 油管 壁面 沥青 
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质 厚 度 的 结论 一 致 ,沉积 能 力 随 着 流速 的 增 大 而 
减 小 。 
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a 


mn 
So 


0.5 1.0 ES 2.0 2:5 
流速 /(m * s?) 
到 12 流速 对 沉积 速率 的 影响 
Fig.12 The effect of oil velocity on deposition rate 


3 结 论 


本 研究 以 水 平井 水 平井 简 内 的 液 固 两 相 流 的 基 
本 理论 为 基础 , 消 流 模型 采用 Realizable k-e 模型 , 固 
体 颗 粒 的 追踪 采用 DPM 离散 相模 块 ,沥青 质 聚 集 颗 
粒 沿 入口 垂直 面 均匀 注入 ,速度 与 油 流 的 速度 保持 
一 致 。 基 于 流 固 单 向 相 耦 合 的 方法 ,模拟 了 流速 
0.3 ~2. 0 m/s 时 不 同 尺寸 沥青 质 颗粒 在 油管 中 的 沉 
Wu. fe. 

1) 在 层 流 阶 段 , 沉 积 百分数 和 速率 随 着 黏度 的 
增 大 而 减 小 ,这 是 由 于 流体 的 流动 平稳 , 管 壁 附近 的 
流体 梯度 低 ,Saffman 升力 减 小 引起 的 ; 滑 流 时 , 当 原 
ith BE 0.001 2 Pa. s 增加 至 0.012 Pa * s 时 , 沉 
积 百 分 数 和 沉积 速率 会 随 之 增 大 ; 2 Jt ut ih RE 
0.012 Pa - s 继续 增加 至 0. 053 Pa - s 时 , Re 从 
10 103 减 小 至 2 300 ,致使 油管 中 的 原油 惯性 力 减 小 ， 
因此 ,沉积 速率 和 沉积 百分数 会 随 之 下 降 。 

2) 在 流速 恒定 时 ,沥青 质 颗粒 沉积 百分数 和 沉 
积 速 率 随 着 颗粒 直径 的 增 大 而 增 大 ,同时 通过 分 析 
出 口 粒 径 颗粒 分 布 得 出 直径 大 于 80 um 的 沥青 质 颗 
粒 大 部 分 都 沉积 在 了 管 壁 ,表明 颗粒 尺寸 越 大 ,沉积 

3) 当 颗 粒 直径 恒定 ,流速 从 0.3 m/s 增加 至 
2 m/s ,沉积 概率 大 概 降 低 了 0.32% , 沉积 速率 减 小 
了 8% 。 这 是 因为 随 着 流速 .阻力 的 增加 ,壁面 产生 
较 大 的 剪 切 应 力 ,使 得 壁面 颗粒 的 苛 附 概率 降低 , 同 
时 也 会 剥离 沉积 在 管 壁 的 沥青 质 沉积 。 
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